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A  general  formalism  is  presented  to  study  hot  carrier  relaxation  in  the  valence  band  of  strained  and 
unstrained  semiconductors.  The  approach  is  based  on  a  six-band  k-p  description  of  the  valence  band 
and  an  anisotropic  Monte  Carlo  method.  We  show  the  results  of  our  studies  for  the  Si-Ge  system. 
Carriers  are  initially  injected  in  the  split-off  band  and  the  carrier  distribution  is  followed  in  time. 
Results  are  presented  for  energy-dependent  energy  relaxation  time.  The  relaxation  times  for  Si  are 
about  10" 13  s  while  those  for  Ge  are  an  order  of  magnitude  higher.  To  study  the  effect  of  biaxial 
strain  produced  through  epitaxy  on  energy  relaxation  times  we  present  results  for  Si0.8Ge0  2  on  a 
{100}  silicon  substrate.  This  work  is  of  relevance  to  the  interpretation  of  pump-probe  experiments, 
which  are  primarily  dependent  on  the  relaxation  of  the  average  of  the  carrier  energy.  ©  1996 
American  Institute  of  Physics.  [S002 1-8979(95)00122-9] 


Hot  carrier  transport  has  been  an  important  issue  in  both 
electronic  and  optoelectronic  devices  and  energy  relaxation 
is  one  of  the  most  significant  factors  in  hot  carrier  transport. 
Theoretical  and  experimental  studies  of  hot  carrier  relaxation 
have  primarily  focussed  on  electrons  in  direct  gap  materials 
such  as  GaAs.1"6  These  studies  have  been  very  valuable  in 
understanding  pump-probe  studies  using  ultrafast  lasers.  It 
has  been  widely  assumed  that  hot  hole  relaxation  times  are 
extremely  fast  and  do  not  affect  the  interpretation  of  the 
pump-probe  studies.  With  continuous  improvements  in  the 
temporal  resolution  of  the  pump-probe  studies,  it  is  becom¬ 
ing  possible  to  study  phenomena  at  the  time  scales  of  — 20 
fs — times  which  are  relevant  to  hot  hole  relaxation.  In  this 
communication  we  present  results  on  hot  hole  relaxation  for 
the  Si-Ge  system.  Not  only  is  this  an  important  material 
system  for  technology,  but  it  provides  us  a  set  of  valence- 
band  structures  where  the  densities  of  states  vary  over  almost 
an  order  of  magnitude,  especially  near  the  band  edges.  The 
system  also  provides  us  an  opportunity  to  study  the  effect  of 
strain  produced  by  epitaxial  growth  on  carrier  relaxation. 

In  modeling  the  valence-band  structure,  we  use  a  six- 
band  k  •  p  Hamiltonian  in  which  the  effect  of  the  strain  can 
be  included  via  the  deformation  potential  theory.7,8  The  ef¬ 
fect  of  alloying  is  included  through  the  virtual  crystal  ap¬ 
proximation  by  a  simple  weighted  average  of  the  Kohn- 
Luttinger  parameters.  This  band  structure  adequately 
describes  the  heavy-hole  (HH),  light-hole  (LH),  and  split-off 
(SO)  bands.  When  Si-Ge  is  grown  on  a  Si  substrate  [along 
the  (001)  direction]  the  strain  removes  the  HH-LH  degen¬ 
eracy  that  occurs  at  k  =  0.  This  is  known  to  reduce  the  near¬ 
band-edge  density  of  states  in  the  strained  alloy.9  In  this 
communication  we  focus  on  the  following  material  systems: 

(i)  silicon; 

(ii)  unstrained  Si0#8Ge0-2; 

(iii)  strained  Si0.8Ge0  2  on  a  {100}  silicon  substrate; 

(iv)  unstrained  Si0i2Ge0  8;  and 

(v)  germanium. 

The  SiGe  system  provides  us  sets  of  valence  bands  with 
the  following  interesting  features. 


(i)  The  valence-band  density  of  states  in  Ge  is  about  an 
order  of  magnitude  lower  than  that  in  Si. 

(ii)  The  strained  alloy  Si0.8Ge02  has  a  sharp  reduction  in 
the  density  of  states  near  the  band  edge  when  com¬ 
pared  with  the  unstrained  alloy;  however,  as  energy  is 
above  0. 1  eV,  the  density  of  the  states  becomes  simi¬ 
lar  to  that  in  the  unstrained  alloy. 

(iii)  Si  has  a  relatively  constant  fractional  density  of  states 
in  the  HH  and  LH  bands  as  a  function  of  the  energy, 
while  Ge  does  not. 

The  consequences  of  these  effects  are  important  to  under¬ 
stand  the  hot  hole  relaxation  in  the  SiGe  system  and  are 
discussed  below. 

The  Monte  Carlo  program  used  for  this  study  includes 
the  effect  of  band-structure  anisotropies  on  the  angular  de¬ 
pendence  of  the  scattering  rates.  The  details  of  the  Monte 
Carlo  method  are  given  in  Refs.  10  and  11.  Monte  Carlo 
simulations  are  performed  using  acoustic-  and  optical- 
phonon  scattering  mechanisms.  Alloy  scattering  is  included 
as  an  additional  scattering  mechanism  in  the  SiGe  alloys. 
Among  these  scattering  processes,  optical-phonon  emission 
and  absorption  play  important  roles  in  energy  relaxation,  be¬ 
cause  the  other  scattering  processes  are  elastic.  The  holes  in 
the  calculation  start  with  the  same  energy  in  the  SO  band. 
Results  are  reported  for  zero  electric  field  and  at  r=300  K. 
During  the  calculation,  average  energy  and  band  occupancy 
are  recorded  as  a  function  of  time. 

In  Fig.  1,  the  band  occupancies  in  silicon  and  germa¬ 
nium  are  shown  as  a  function  of  time.  Carriers  start  in  the  SO 
band  with  an  initial  energy  of  0.5  eV  for  both  material  sys¬ 
tems.  Two  observations  must  be  made  from  the  results 
shown: 

(i)  The  carriers  transfer  out  of  the  SO  band  very  rapidly 
in  both  Si  and  Ge; 

(ii)  for  Si  the  ratio  of  the  occupancy  between  the  HH  and 
LH  bands  reaches  a  constant  value  long  before  the 
carrier  reaches  equilibrium  energy  (this  result  is  dis¬ 
cussed  in  Fig.  2,  below). 
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FIG.  1 .  Band  occupancy  of  each  band  in  silicon  and  germanium  as  a  func¬ 
tion  of  time.  The  results  shown  correspond  to  cases  where  holes  with  energy 
of  0.5  eV  are  placed  in  the  split-off  band. 

On  the  other  hand  for  Ge,  the  HH  and  LH  band  occupation 
reaches  a  constant  value  only  when  the  holes  have  reached 
their  equilibrium  energy. 

To  understand  these  observations  we  examine  some  of 
the  scattering  rates  important  for  carrier  relaxation  and  en¬ 
ergy  loss.  We  find  that  holes  in  the  SO  band  of  Si  transfer  to 
HH  with  an  interband  scattering  rate  that  is  2.5  times  larger 
than  that  in  Ge,  due  to  the  high  density  of  states  in  the  silicon 
HH  band.  At  0.5  eV,  the  interband  scattering  rate  from  SO  to 
HH  in  Si  is  approximately  1X1014  s_1  and  is  an  order  of 
magnitude  greater  than  the  SO  to  LH  scattering  rate,  and  two 
orders  of  magnitude  greater  than  the  SO  intraband  scattering 
rate.  Because  of  the  high  rates  of  interband  transfer,  holes  are 
transferred  to  the  HH  band  in  approximately  50  fs  in  Si.  In 
Ge  this  time  is  a  little  longer  but  it  is  still  less  than  100  fs. 

To  understand  the  second  observation  noted  above,  we 
examine  the  ratio  of  the  density  of  states  (DOS)  of  the  HH 
and  LH  bands  as  a  function  of  energy  in  Ge  and  Si.  Figure  2 
shows  the  ratio  of  DOS  in  each  band  to  the  total  DOS.  At 


FIG.  2.  Ratio  of  the  density  of  states  in  each  band  to  the  total  values  in 
silicon  and  germanium  as  a  function  of  energy. 


FIG.  3.  Average  energy  in  silicon,  unstrained  and  strained  Si08Ge02,  un¬ 
strained  Si02Ge08,  and  germanium,  (u)  and  (s)  represent  unstrained  and 
strained  cases,  respectively.  In  all  cases,  the  holes  are  initially  placed  with 
0.5  eV  energy  in  the  split-off  band. 


high  energies,  Ge  has  roughly  the  same  ratio  of  DOS  in  HH 
as  in  silicon  but  as  we  move  closer  to  the  band  edge,  while 
the  Si  ratio  remains  energy  independent,  in  Ge  it  increases  at 
0.1  eV  and  converges  to  a  higher  value  than  for  Si.  There¬ 
fore,  the  band  occupancy  in  Si  reaches  equilibrium  value  just 
after  the  transition  from  SO  to  HH,  while  that  in  Ge  takes  a 
much  longer  time  and  comes  to  equilibrium  only  when  the 
hole  energy  has  thermalized. 

In  Fig.  3  the  average  energies  in  Si,  unstrained  Si0  2Ge0  8, 
strained  Si0.2Ge08,  unstrained  Si08Ge02,  and  Ge  are  pre¬ 
sented  as  a  function  of  time.  The  initial  energy  is  0.5  eV  for 
all  cases.  Due  to  the  large  scattering  rates,  holes  in  Si  lose 
their  energy  much  faster  than  those  in  Ge.  The  unstrained 
SiGe  alloys  have  intermediate  values  between  them,  since 
energy  relaxation  is  closely  related  to  the  DOS.  The  holes  in 
the  strained  Si02Ge08  initially  lose  their  energy  with  the 
same  rate  as  for  the  unstrained  system;  however,  they  have 
larger  relaxation  time  as  their  energy  decreases.  This  is  be¬ 
cause  the  strained  system  has  a  lower  effective  mass  at  band 
edge  due  to  the  strain  induced  band  deformation;  therefore, 
the  effect  of  the  strain  can  be  seen  when  the  average  energy 
is  low.  Results  for  other  SiGe  compositions  can  be  inferred 
from  the  cases  discussed  in  this  communication  since  the 
relaxation  times  are  closely  related  to  the  density  of  states  of 
the  alloy. 

The  results  shown  in  Figs.  1  and  2  are  for  carriers  ini¬ 
tially  injected  at  0.5  eV.  As  the  carriers  relax,  the  energy 
distribution  is  no  longer  sharp  but  starts  to  spread  around  an 
average  energy.  To  estimate  the  effect  of  this  energy  spread, 
we  calculate  an  average  energy  dependent  relaxation  time. 
We  define  this  relaxation  time,  r,  through  the  temporal  gra¬ 
dient  of  the  average  energy  ( E ), 

1  _  1  d(E ) 

T«£))-~(£)~rfT '  (1) 

In  Fig.  4  the  dependence  of  the  energy-dependent  relaxation 
time  on  the  average  energy  is  given  in  silicon  and  germa- 
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FIG.  4.  Dependence  of  the  energy  relaxation  time  on  the  average  energy  in 
silicon  and  germanium.  The  relaxation  time  is  defined  by  the  temporal  gra¬ 
dient  of  the  average  energy.  The  results  are  shown  for  three  different  initial 
energies. 

nium.  We  use  three  initial  conditions  for  these  results.  Holes 
are  injected  into  the  SO  band  with  initial  energies  of  0.3,  0.4, 
and  0.5  eV.  We  see  that  even  though  the  initial  conditions  are 
different  in  the  three  cases,  the  rate  of  energy  relaxation  is 
dependent  on  the  average  energy  and  is  not  significantly  af¬ 
fected  by  the  spread  in  the  energy  distribution. 

In  conclusion,  energy  relaxation  in  the  valence  band  has 
been  calculated  using  Monte  Carlo  simulations  with  a  k  •  p 
band  structure.  We  have  shown  results  on  carrier  transfer 


times  from  the  SO  band  to  the  HH  and  LH  bands  as  well  as 
the  time  evolution  of  the  band  occupancy.  The  effect  of  den¬ 
sity  of  states  and  the  influence  of  strain  on  carrier  relaxation 
have  also  been  discussed.  We  expect  that  the  results  pre¬ 
sented  here  will  be  useful  for  understanding  ultrafast  laser 
pump-probe  measurements  focusing  on  hot  hole  relaxation. 
In  such  experiments  it  is  the  average  carrier  energy  that  is  of 
importance,  more  than  the  details  contained  in  the  shape  of 
the  distribution  function.  Our  results  show  that  the  intrinsic 
properties  of  the  relaxation  rates  are  such  that  its  value  at  a 
given  energy  is  virtually  independent  of  the  hole’s  initial 
energy. 

This  work  was  supported  under  the  Army  URI  program 
(DAAL  03-92-G-0109)  and  under  AFOSR  Contract  No. 
AFOSR  91-0349. 
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